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We present a thorough description of the physical regimes for ultracold bosons in double wells, 
with special attention paid to macroscopic superpositions (MSs). We use a generalization of the 
Lipkin-Meshkov-Glick Hamiltonian of up to eight single particle modes to study these MSs, solving 
the Hamiltonian with a combination of numerical exact diagonalization and high-order perturbation 
theory. The MS is between left and right potential wells; the extreme case with all atoms simulta- 
neously located in both wells and in only two modes is the famous NOON state, but our approach 
encompasses much more general MSs. Use of more single particle modes brings dimensionality into 
the problem, allows us to set hard limits on the use of the original two-mode LMG model commonly 
treated in the literature, and also introduces a new mixed Josephson-Fock regime. Higher modes 
introduce angular degrees of freedom and MS states with different angular properties. 



I. INTRODUCTION 

The superposition of quantum states Is a postulate in 
quantum mechanics. When these states can be distin- 
guished macroscopically it leads to the fundamental ques- 
tion of how the theory describing the physics of point par- 
ticles or atoms is connected with macroscopic objects [1]. 
Ultracold bosons in double wells provide a useful realiza- 
tion of macroscopic superposition (MS) states, in which 
the distinguishable macroscopic property is the localiza- 
tion of atoms in one of the wells [2]-(4]. In this system, it 
is possible to restrict the atoms to occupy only two single 
particle modes, corresponding to their being condensed 
in the lowest energy state and localized spatially in one of 
the wells [5ll^. In this simplified two mode picture, the 
main two processes are the tunneling of atoms between 
wells and their interaction in pairs in one of the wells. 
When the tunneling energy is larger than the interaction 
energy, this is an ideal system for studying Josephson 
effects [71119). Conversely, when interactions dominate 
over tunneling, MSs can be obtained [H 123 HH • When 
all N atoms occupy simultaneously the two single parti- 
cle states localized in each well, these MSs are known as 
NOON states. The tunneling dominated regime is known 
as the Josephson regime, while the interaction regime is 
known as the Fock regime. In the latter, the two mode 
approach can cease to be sufficient to describe the sys- 
tem, due to large interactions populating single particle 
excited states. Here, we offer a detailed study of the 
eigenvector spectra in one dimensional (ID) and three 
dimensional (3D) double wells, with special attention to 
the appearance of MS states which occur in the Fock 
regime, and we include more modes in our approach, as 
appropriate for this regime. 

To model ID double wells we use a four-mode gen- 
eralization of the Lipkin-Meshkov-Glick Hamiltonian 
(LMGH) [551 123] > where the atoms are allowed to occupy 
two excited on-well (sometimes called on-site) localized 
single particle eigenstates (see Fig.[TJa)). Due to the new 
two excited modes, there are three new processes in the 




FIG. 1. (Color online) Schematic of the Double Well, (a) 
ID: the J's indicate single-particle tunneling processes while 
the f/'s indicate two-particle tunneling and interaction pro- 
cesses, (b) 3D: the spherical harmonics Y^™, with total and 
z-component of the angular momentum £ and m respectively, 
are distorted by the coupling to the nearby well. 



model. On the one hand, pairs of atoms in the ground 
state in one well can be excited to the first level, or de- 
cay from this level to the ground state, as sketched in 
Fig. [l](a). Also, the atoms excited to the first level of 
eigenenergies can interact on-well or tunnel to the other 
well, also sketched in Fig. [T](a). Then, in this picture, 
the conventional division into the Josephson and Fock 
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regimes must be amended to include a new mixed regime, 
in which the atoms behave differently depending on their 
energy level, as discussed in Sec. |Vj We use perturbation 
theory to characterize the eigenstates and eigenvalues in 
the Fock regime and we also characterize them analyt- 
ically in the non-interacting regime. In between these 
two regimes, we use numerical exact diagonalization of 
the LMGH to show how the transition between these 
regimes occurs. We identify clearly the regimes in which 
the original NOON states as well as novel MS states can 
be obtained. We characterize the limits of the model, the 
occurrence of crossings in the spectra, and couplings be- 
tween states with atoms in different levels. Our crossing 
criteria clarify when the one-level approximation (LMGH 
with only two modes) ceases to be sufficient, and one 
must use at least the two-level approximation. 

The four-mode, two-level generalization of the LMGH 
is the natural way to introduce dimensionality in the 
problem. When the same two levels are generalized to 3D 
double wells, an eight-mode generalization of the LMGH 
is obtained as sketched in Fig. [l|b). Then, a new quan- 
tum number, the orbital degree of freedom, is introduced 
and it plays a key role [23]. Since new parameters arise 
in the 3D case, new dynamical regimes are found. More- 
over, new MS states involving angular degrees of freedom 
can occur. In a previous work we characterized the lim- 
its of the one- and two- level approximations, the occur- 
rence of crossings in different regimes, and the coupling 
between states with atoms in different levels [24] . In this 
Article we review, complete, and extend the ID picture 
obtained in [31 [53]. We also compare with the 3D case, 
and summarize the derivation of the 3D two-level LMGH 
and the routes to new physics in 3D which are currently 
being explored. 

By means of a semiclassical, mean-field approximation 
macroscopic quantum tunneling has been predicted in 
the Josephson regime, while macroscopic self-trapping in 
one of the wells has been shown to occur when the in- 
teractions grow [7l416|. The macroscopic quantum tun- 
neling and self-trapping were observed in a recent ex- 
periment [17] . The former can be understood as the 
d.c. Josephson effect for superconductors, while the a.c. 
Josephson effect corresponds to small oscillations in the 
difference between the populations in both wells obtained 
for bigger interactions [TH]. The a.c. Josephson effect was 
recently demonstrated experimentally [19] . One perspec- 
tive is that the energy levels depend on the occupation 
of one of the wells. Then one can analyze the double well 
problem in terms of Landau-Zener tunneling. In this case 
the energy between levels changes with a constant small 
rate, leading to an adiabatically driven Josephson junc- 
tion [26ng. 

In the Fock regime the semiclassical approach ceases 
to be appropriate, and other approaches can be useful; 
the most numerically accurate but also the most compu- 
tationally intensive is the multiconfigurational Hartrce 
method [29 33 . Methods based on the LMGH [221 [25] . 
although limited to weaker interactions, have the advan- 



tage that one can treat more particles in higher dimen- 
sions, as well as make use of perturbation theory to ob- 
tain analytical results [21 [M] . In this Article we use a 
modified LMGH for these reasons. The LMGH was first 
introduced in the framework of nuclear physics [22] . It 
has applications in a wide range of other fields, like ul- 
tracold bosons, and it is a perfect simple Hamiltonian 
to study a vast range of quantum effects, from quantum 
phase transitions [23 EH] to Josephson oscillations. The 
LMGH also applies in cold atoms for two hyperfine boson 
species in a single well, where MS states are predicted to 
occur [321 [37]. In this Fock regime the two mode ap- 
proach could be insufficient to characterize the problem, 
since crossings and coupling to other more excited single 
particle states may be possible. Indeed, it is necessary 
to consider a second level in ID for the physical param- 
eters found in typical BEC experimental systems [2] . To 
consider this possibility we use a four-mode generaliza- 
tion of the LMGH for ID double wells and an eight mode 
characterization in the 3D case. 

Besides being a good system for fundamental studies 
of quantum many body effects, ultracold bosons in dou- 
ble wells also have many technological applications, for 
example, in quantum high precision measurements of in- 
ertial and gravitational fields [351 ISH] , as a primary ther- 
mometer [40] . or in quantum computing f 41H43j . 

The Article is organized as follows. In Sec.[II|we intro- 
duce the initial second quantized Hamiltonian. In Sec.|III| 
we obtain the ID four-mode LMGH. Then, in Sec. |IV| we 
characterize the eigenvectors and eigenvalues in ID for 
the Josephson and Fock regimes. In Sec. |V| we obtain 
the limits of the model and the occurrence of crossings, 
we discuss all possible regimes, and we offer numerical ex- 
amples to illustrate the results. Finally, we derive the 3D 
Hamiltonian and compare with the ID case in Sec. |VI[ 
We conclude, summarize, and discuss lines of future re- 
search in Sec. IVIII 



II. SECOND QUANTIZED HAMILTONIAN 

The second quantized Hamiltonian for N interacting 
bosons of mass A4 confined by an external potential ^(x) 
with X S M'^ in terms of the bosonic creation and anni- 
hilation field operators ^^(x), ^^^(x), obeying the usual 
commutation relations, is 



H = y'dx*t(x) _-^v2 + V^(x) *(x) 
+ i ydx4't(x) J dx'*t(x')X^„t(x-x')'I'(x') 



(1) 



where t^nt(x — x') stands for the two-body interaction 
and V{x) is a three dimensional double- well potential 
with minima at x = ±a G and a local maximum at 
x = 0. Without loss of generality, we assume V{x) = 
V{x) -\- V{y) -I- V{z). We take the ID potentials in the 
X and y directions as isotropic harmonic oscillators, that 
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is, Vx 



^ujIx"^ and Vy 



2^x-^ •'y — h^lv^ ^ where both frequencies 

are of the same order. For simphcity, we consider lOx = 
Ljy = u)_\_ in the fohowing. A conventional ID double well 
potential V{z) of barrier height Vq can be approximated 
near its minima at z = ±a by 



V{z±a) w ]^Ml^^z^ . 



where 



1 d'^V 
M dz^ 



1/2 



(2) 



(3) 



is the local trapping frequency in each well. The recoil 
energy, defined as Er = h? /Ma?, is used through the 
dimensionless parameter V^/Er to determine the barrier 
height in most of the experiments, with A = 2a an effec- 
tive wavelength. 

We assume first that uj is of the same order as uj±^ and 
we consider low densities, for which only binary colli- 
sions are relevant; also, we consider low energies, when 
these collisions are characterized by the s-wave scatter- 
ing length of the atoms, a^. The diluteness condition for 
a weakly interacting Bose gas is 



< 1. 



(4) 



where n is the average density of the gas. In the context 
of the double- well potential, the maximum density of the 
gas is approximately n — N / {^/2tt aho)^ : where Oho is the 
harmonic oscillator length given by ajjo = y^h/Aiuj (a 
similar condition must hold in the other two directions). 
Correspondingly, we restrict our discussion to the regime 



7V1/3 ^ 



27r aho/a. 



(5) 



Although the system is said to be weakly interacting 
when condition ([5| is met, the interaction energy can 
be on the order of the kinetic energy and dilute gases 
can therefore exhibit non-ideal behavior [SJ IHj- Under 
these conditions the second quantized Hamiltonian can 
be approximated as 



2M 



V{x) 



*(x) 



d3x#t(x)*^(x)^(x)#(x). 



(6) 



The coupling constant g is proportional to the s-wave 
scattering length, g — Anfi^as/M; in Eq. ^ we took the 
two-body potential from Eq. ([T]) to be approximated by 
an effective local interaction V^nt(x — x') ~ g S^^^ (x — x'). 



III. ONE-DIMENSIONAL 
LIPKIN-MESHKOV-GLICK HAMILTONIAN 

Let us consider first the reduced ID case, taking w <C 
The double- well potential can be reduced to one 



spatial dimension in extremely anisotropic traps, where 
the transverse trapping frequencies must be sufficiently 
high to reduce the dimensionality of the single-particle 
wavefunctions. However, one must avoid potential reso- 
nances by not squeezing the trap too tightly [U]: it is 
sufficient t hat the transverse harmonic oscillator length 
Oho, J- = \/fi/M.LiJj_ 3> ds- Under these assumptions, we 
can reduce the dimensionality of the second quantized 
Hamiltonian, by considering that the particles interact 
and tunnel only in one-dimension. Then, we can use 
appropriate superpositions of the eigenfunctions of the 
one-dimensional single particle Hamiltonian 



(7) 



to obtain a set of on-well localized wavefunctions, in an 
analogous way that on-site Wannier states are obtained 
from Bloch functions on a lattice \'> . We use this set of 
functions to expand the field operator: 



^(z) =^bjeijji{z - zj), 

3,1 



(8) 



where ipe{z—Zj) are the on- well localized functions. Here, 
zl = —a and zr = a are the minima of the left and right 
wells, respectively, j G {L, R} is the well index, and the 
label £ is the level index. The level index increases with 
single particle energies in each well. For the two level 
approach considered here, £ € {0,1}. The operators bje 
and b^jf satisfy the usual bosonic annihilation and cre- 
ation commutation relations. 



[bji, bj, 



0. 



(9) 



Let 'p'^{z) be the nth eigenfunction of Hamiltonian ([7|, 
with eigenvalue e„. Then, the localized functions at well 
j are 



iPjiiz) 



^{^\z)±^^iz)) 
^{^^z)±^\z)) 



(10) 
(11) 



with j — L for the plus sign and j = R for the minus sign. 
The corresponding eigenvalues are Eq ^ [e^ + e^) / 2 and 
El — {e^ + e'*)/2. These on-well localized eigenfunctions 
are represented schematically in Fig. [l](a) . 

Substituting Eq. ([s]) into the second quantized Hamil- 
tonian ([b]) yields the two-level Hamiltonian: 



H — Hq + Hi + Hq 



(12) 



with 



Hi = -Ji ^ b^'^ibj'i + Uu ^ fiji {fiji - 1) + E'^ ^ fiji , 

(13) 
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(14) 



where fije = b^j^bjg is the number operator. The hopping 
and interaction terms are 



Ji 



and 



dz^pl (z — a) 



2M 



Viz) 



Mz+a): (15) 



I / dz\Mz)\'\i:Az)\' 



(16) 



respectively. The coupUng constant obtained after reduc- 
ing the dimensionaUty of the problem is gi = 2huj±as. In 
the following, we denote Uu simply as U^. 

The functions ipii^) resemble roughly the eigenfunc- 
tions of the harmonic oscillator potential: 



1/2 -1/4 



V'o(z) = Oho 



V2aL 



(17) 



V'i(z) = ai/'V27r-i/4(zKo)e-"'/'"'" . (18) 

The energy of an atom associated to the wavefunction 
i\)t(z — Zj) is El K. T(u>(l + 1/2). This approximation 
fails when the overlap between different on-well localized 
eigenf unctions is small. Then, in general, it gives a poor 
approximation of the Jg coefficients given in (15). Never- 



theless, they can be used to obtain a good approximation 
for the overlap integrals purely on-well, like those defin- 
ing the coefhcients, Eq. (16). This is valid for both 



regimes, though it gives only a rough idea of these coeffi- 
cients when the actual single particle eigenfunctions are 
highly distorted, as occurs in the Josephson regime. In 
this case, it can only be used to obtain scaling relations 
in the problem. By solving these integrals analytically 
using Eq. ( 18 ) to approximate the single particle eigen- 



functions, we obtain 



gi /aho\ 
J2^ \ a J 



(19) 



Moreover, we can relate the other interaction coefhcients 
to Uq by solving the corresponding integrals, which gives 
Ui = (l/2)C/o and Uoi — (3/4)[/o- Here we use these ex- 
pressions in the illustrative examples presented in Sec.jVj 
though for actual potentials to obtain accurate values 
for these coefficients, it is necessary to calculate the sin- 
gle particle eigenfunctions numerically to solve the inte- 
grals (16). 



Finally, to deduce the Hamiltonian ( 12 ) we neglected 
off-site interactions. These correspond to interaction 
terms like: 



TI" - 



51 
2 



dz\il}fXz - Zj)f\'4!e,{z - zy)Y 



(20) 



where j ^ j'. If the barrier is infinitely high 
these integrals vanish. If the single particle eigenfunc- 
tions are approximated by Eq. ( 18 ) this integral gives 



?7o exp[— (a/aho)^]- Then, one could assume that these 
terms can be neglected for sufficiently high or separated 
wells, that is, if Oho ^ a- Nevertheless, these are not on- 
well integrals and then the expressions of the harmonic 
oscillator eigenfunctions, Eq. (18 1 do not give accurate 



results. Then these integrals have to be evaluated numer- 
ically for the actual potential and for the correct numer- 
ically evaluated single particle eigenfunctions. Moreover, 
these coefficients have to be compared with Ji to justify 
their being neglected, since both quantities could be com- 
parable in the Fock regime. In the Josephson regime we 
can safely neglect these coefficients since they are much 
smaller than any other interaction coefficient, which in 
turn are smaller than the tunneling ones. For the Fock 
regime, we assume that the tunneling coefficients are al- 
ways bigger that these cross-terms, and have verified that 
our assumption holds for a number of specific cases. We 
refer the reader to Ref. [46^ for a thorough discussion of 
the effect of these terms. 



IV. CHARACTERIZATION OF THE 
EIGENSTATES 

A. General Considerations 

An arbitrary state vector |^), in the solution space of 
the two-level Hamiltonian Eq. (12), can be expressed in 
Fock space as 



where 



n-i 

i=0 



|0) 



(21) 



(22) 



Here is the dimension of the Hilbert space i 
is the Fock index, and |cip is the probability of finding 
the atoms distributed between different energy levels and 
wells in the Fock vector when the system is described 
by state |^'), that is q = (*|^); henceforth we refer to 
Ci as the Fock- space amplitude. The Fock vector \i) can 
also be written in longer form as 



LO' "'flO/l"'Ll' "-Rll ' 



(23) 



indicating the occupation of each well (L,R) and each 
level (0,1) associated with Fock index i. We work in the 
canonical ensemble, i.e., we require the total number of 
particles 



(24) 



to be constant. 

The Fock index i is chosen to increase starting with the 
number of atoms in well j = L, and then subsequently 
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increasing with the number of atoms in the same, left 
well moving up into the excited level i = 1. Therefore, 
for the first N + 1 Fock vectors i — 1 + rij^Q and they 
correspond with vectors with no occupation of the excited 
level. Then, they satisfy 



1 



(St.)' 



'm 



(25) 



for i = 0, 1, • • • , N . For example, for iV = 2, the Fock 
vectors with index i = 1,2,3 are |0,2)|0,0), |1,1)|0,0), 
and |2, 0)1 0,0), respectively. The one- level approxima- 



tion can easily be recovered from Eq. (12 1 by requiring 
i < N + \. In this truncated space, the dimension of the 
Hilbert space reduces to that of the one-level approxima- 
tion, namely, iV -I- 1, and the two-level Hamiltonian H 
reduces to the one-level Hamiltonian Hq. 
A general expression for the Fock index is: 

Afi-l 

i = l + nLo+ {N +l-p){l+p) 

+ {N + 1- Ni)nLi, (26) 

where Ni is the number of atoms in the excited level. 
Then, the following Fock vectors for iV = 2 are |0, 1) |0, 1) , 
|1,0)|0,1), |0,1)|1,0), and |1,0)|1,0) with indices i = 
4, 5, 6, 7. The last three show occupation only of the 
excited level and they are |0,0)|0,2), |0,0)|1,1), and 
|0,0)|2,0), for i = 8,9,10. Thus for just two atoms 
there are already 10 states in our two-level problem, and 
with increasing A'^ it is necessary to use numerical matrix 
methods to keep track. 

The eigenstates l^^*'-') of the two- level Hamiltonian 



(121 satisfy: 



(27) 



where e^'^-' is the energy eigenvalue corresponding to the 
state |(/>''^-'). The eigenstate label k is chosen to increase 
with e'^'^). In order to describe these states, we will use 
the previously introduced Fock-space amplitudes 



(28) 



now containing two separate indices for clarity: i is the 
Fock index, describing the ordering of the Fock basis; k 
is the energy index, describing the ordering of the energy 
eigenvalues. 

When interlevel effects are not relevant, the eigenstates 
fall into one of two categories: harmonic-oscillator-like 
states, and macroscopic superposition states (MS states). 
When the barrier between wells is low, Jq ^ NUq, all 
states are harmonic oscillator-like (HO) states. This is 
the Josephson regime and it is characterized by 



U = N/Ci « 1 , 



(29) 



where Q = Ji/Ui. This criterion must be evaluated sepa- 
rately for both levels, since Ji > Jq, while the interaction 
terms are of the same order. 



On the other hand, MS states dominate the spectrum 
in the high barrier limit. This is the Fock regime, in 
which 



(30) 



Again, this criterion must be evaluated separately for the 
two levels. 

When the level spacing is comparable to NUq, inter- 
level effects can no longer be neglected, and a third cat- 
egory of eigenstates emerges. These states show weak 
coupling between states with particles only in the lowest 
energy level and others with particles in the excited one. 
We name these coupled excited states shadows of the MS 
states, because in a surface plot of the Fock-space ampli- 
tudes as a function of both energy and Fock index, they 
appear as faint copies of the MS states with occupation 
only of the lower level at higher Fock index [24j. In the 
following two subsections we characterize the eigenstates 
in the Josephson and Fock regimes. 



B. Non-interacting Regime 

Let us consider the extreme Josephson regime, or non- 
interacting regime, for which Ui = 0. The Hamiltonian 
reduces to H = Hq + Hi, with 

Hi = -Ji ^T^r +^iY. ■ (31) 

Since the Hamiltonian is clearly separable, its eigenstates 
l^*^*^)) are a direct product of the one-level ones, j^^^^f '>): 

|0('^))=(g)|0(^^''')). (32) 



where K^'^^ = 0, 1, . . . , N^'^' is the one- level eigenstate la- 

(k) 

bel and for the fcth eigenstate there are A'^ atoms at 
level £. Since the total number of atoms is conserved, 

)ne-level eiger 

(fe), (fe) Jk). 



pressed as {(p^^i '') — J^^it'l''^^ ^''^Rf^ with 



J2i=o 1 ^e' ~ one-level eigenstates can be ex- 



(fc)^ 



(k)_ 



0K7 



A^, 



(fe) 



(33) 



where n^*^^ is the occupation of the left well for the corre- 



sponding Fock vector, h^{k}{n^^^\N^''^) is a A'^'^^'th order 

is the square 



,(fe). 



discrete Hermite polynomial, p{n^^^ 
root of the binomial distribution, 

1 



A, 



(fe) 



p (i 

and Oj-ik) {Ng'^'' 



2wf' /2 



\ 



Ak) 



ni,!(A; 
is the normalization factor 



(34) 



a^(.)(A^f)) 



(fe) 



if, 



ik)^ 



(35) 
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The eigenenergies e^''^ can be expressed in terms of the 
one-level eigenenergies as 



where 



--ik).i 



= -Jf{Nl''^-2K. 



(k) 



(36) 



(37) 



The one-level Hamiltonians, expressed in the Fock ba- 
sis is, in the non-interacting limit, a harmonic oscillator 
potential truncated at hard walls. The analytical ex- 

(k) 

pressions of the Fock-space amplitudes c-^ resemble the 
probability amplitudes of observing the particle in posi- 
tion X for the ID harmonic oscillator problem, obtained 
from the expression of its eigenfunctions in the position 
representation. Moreover, in both cases the eigenvalues 
Eke are linear in K^'^K The analogy can be established 
between occupation of the wells in the double well and 
the positions in the Harmonic oscillator problem. Since 
one variable is discrete while the other is continuous the 
analogy is valid for the limit of infinite particles. Let us 
show this for the one level approach, which w e ob tain 
considering H = Hq, where Hq is given by Eq. (31). In 



the Heisenberg picture the evolution of the operators bjo 
for j — L, R is given by 



db 



dt 



'j'O ■ 



(38) 



with j' ^ j, and where we considered Eq = for simplic- 
ity. In the following we omit the level index £ because 
we are considering only the one level approximation. For 
N ^ oo, we can approach the operators by c- numbers 
aj = y/Nje^^^. Substituting we obtain the two corre- 
sponding equations of motion which in turn can be ob- 
tained from the pendulum Hamiltonian 

n = -2Jy/NLNRC0s{eL - Or) = -JNy/l - cos{9) , 

(39) 

where iV = Nl+Nr, z = (Nl-Nr)IN, 9 = Ol^Or. The 
previous approach is widely known to be extended to the 
small interacting regime, where it gives a Hamiltonian 
analogous to that of a non-rigid pendulum. The study 
of the dynamics of such a Hamiltonian permits one to 
predict macroscopic quantum tunneling and self-trapping 
within the semiclassical approach [3{9]. In our case, this 
can be extended to a harmonic oscillator Hamiltonian, 
where the angular variables can take any value. The 
analogy here is with the harmonic oscillator truncated 
at hard walls, because the number of particles is finite. 
Due to this analogy, the eigenstates (32) are said to be 



harmonic-oscillator-like. In Sec.|V]we obtain numerically 
the eigenvectors and eigenvalues for small interactions, 
showing that for non-zero interactions the eigenvectors 
and eigenvalues closely resemble the ones characterized 
for the non-interacting regime. 



C. Fock Regime 

In the extreme Fock regime, or the infinite-barrier 
limit, Jo = Ji = 0. In this regime, the coefficient J7oi 
cannot be neglected. Then, the eigenvectors of Hamil- 



tonian (12 1 are Fock vectors, if we further assume that 
2/SE ^^>TPUq. The latter assumption is needed because 
the term bjibjibji'bji' couples Fock vectors with atoms in 
different levels. As demonstrated in Sec.|Vj this coupling 
is small under that assumption. In this limit, the eigen- 
values are 



N, 



(k) 



(k) 

I 

{k) 



2 n 



[k) 



(k)' 



{k) (fe) 



(40) 



Here, the spectra reflects the main symmetry in 
the problem: that of discrete rotations of an angle 
TT around an axis perpendicular to the x — y plane 
that intersects this plane in the origin. This in turn 
means that the eigenstates, which are the Fock vectors 
\nLo,nRQ)\nLi,nRi), are degenerate in pairs with those 
obtained for L ^ R. Small but non-negligible hopping 
terms, Jq <C Uq and Ji <C Uq, break this degeneracy, as 
can be shown by non-degenerate high-order perturbation 
theory. Then, in the small tunneling limit, that is, in the 
Fock regime, the eigenvectors are quasi-degenerate sym- 
metric and antisymmetric superpositions of those pairs 
of Fock vectors that are degenerate in the infinite barrier 
case. For the particular cases in which all atoms occupy 
the same level the eigenstates are direct products of one- 
level vectors 10^)10, 0)^' with 



± 



V2 
1 



1 



'Li 



'Re 



|0) 



N- 



(41) 



for < < N/2. These are MS states in which i' and 
N — V atoms simultaneously occupy the £th energy level 
of both wells; although the expression might appear com- 
plicated, in fact it is just a two-state approximation. We 
have neglected terms on the order of ( J^/C/^)^"^" and 
smaller. Here is the usual arbitrary phase associated 
with vectors in a Hilbert space. We will set 9q — Q for 
the rest of this Article. 

The special case v = represents an extreme MS state, 
or NOON state, in which all N atoms simultaneously oc- 
cupy the left and right wells, which can be written also as 
(1/^2) (|iV,0)|0,0) ± |0,iV)|0,0)). These symmetric {+) 
or antisymmetric (— ) MS states are nearly degenerate, 
with a splitting Ae^(i^) between the states 10^) and 
given by 
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(42) 



up to {N — 2i/)th order in Ji/Ui. General MS states are 
symmetric/antisymmetric superpositions of Fock vectors 
with atoms in both levels, that is, \v,N() — v)\v'^Nx — 
v') and jA^o ^ v,v)\Ni — v\v'), with splittings that are 
proportional to ( Jq/C/q)^""^'' {Ji/Ui)^^'^"' . 

Since Ji > Jq, it is also possible that Jq <C Uq but 
Ji > Uq; this is the mixed regime. The atoms in the 
lower level behave as in the Josephson regime while the 
ones in the excited level behave as in the Fock regime. 
We will show numerical examples of the mixed regime in 
the following section. 



V. CHARACTERIZATION OF THE BOUNDS 
OF THE MODEL, REGIMES, AND CROSSINGS 

A. Limits of the Two-level Model 

Let us find criteria in terms of the characteristic pa- 
rameters of the problem to describe the limits of all 
regimes, the occurrence of shadows of the MS states, 
crossings, and the limits of applicability of the Hamil- 



tonian, Eq. ( 12 ). We begin by seeking criteria of validity 
for the model. On the one hand, the consideration of 
quantum tunneling associated with a second level of en- 
ergies is meaningful if Vq > Ei, as evident from Fig. [TJa) . 
If the eigenfunctions are approximated by the eigenfunc- 
tions of the harmonic oscillator, then this criterion can 
be approximated by Vq > {3/2)huj. Secondly, the dilute- 
ness condition given by Eq. ([s]) can be written in three 
dimensions as N^^^U^^/AE < 1 [24]. In ID, this be- 
comes 



LUj_AE 



< 1 



(43) 



which, contrary to the 3D case, accounts also for the ge- 
ometrical compression in one of the dimensions through 
the term uj/uj±. Here, it is also important to consider the 



restrictions on ujj_ discussed at the beginning of Sec. Ill 
Finally, a last limit of the model is obtained when the 
coupling between levels is big enough to require more 
levels to characterize the eigenvectors. We obtain below 
the criterion delimiting when this coupling is too big. 
If this criterion is not fulfilled, more levels are required 
in the approximation or one should use the MCTDH 
method pSHST] . depending on numerical efficiency re- 
quirements and spatial dimensionality. Other methods 
such as matrix product state algorithms, for example, 
time-evolving block decimation, can also of utility in this 
regime [4 7) . 
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FIG. 2. (Color online) Josephson regime example, (a) Fock- 
space amplitudes for all the eigenvectors, which show the 
harmonic oscillator behavior, (b) Zoom of panel (a), (c) 
Fock-space amplitudes for the ground state, (A'^ -I- l)th, and 
{N + 2)th excited eigenstates. ((d) Eigenvalues as a function 
of Fock state index. Inset: lowest level approximation, first 
N+1 states only. The energy increases linearly, with charac- 
teristic jumps occurring each time a particle enters the excited 
level. 



8 



N'=10 
N =9 




(a) 



50 
40 
30 
20 
10 



50 100 150 

(b) 

0.2 



200 250 ^ 

(c) 



N^=1 
N^=0 





FIG. 3. (Color online) Intermediate regime example. Panel 
layout as in Fig. [2] The lowest excited eigenvectors show 
no occupation of the lower level. Among these, the lowest 
excited 5 show harmonic oscillator behavior, while the next 6 
eigenvectors behave as in the Fock regime. The eigenvalues 
show also this mixed HO-MS behavior. Notice that the A'' + 1 
more excited eigenstates of the spectra are harmonic oscillator 
like. 
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FIG. 4. (Color online) A typical example in the Fock regime. 
Panel layout as in Fig. [2] Here, the eigenstates appear in 
quasidegenerate pairs of MS states. 



B. Discussion of the Different Regimes 

The Josephson and Fock regimes are characterized by 
the criteria discussed above, Eqs. (29) and (30). Next 



we will show some numerical examples of the direct di- 
agonalization of the Hamiltonian ( 12 ) for different sets of 
parameters, corresponding to different regimes. In Fig. [2] 
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Fock index i 


iVo 


iVi 


1 to iV + 1 


N 





TV + 2 to 3Af + 2 


N - 1 


1 


3iV + 3 to 






3iV + 3 + 3(iV - 1) 


TV- 2 


2 



TABLE I. Fock indices for the the Fock vectors up to the 
{?>N + 3 + 3(A'' — l))-th resulting from the ordering given in 
Eq. (26 1. These include all combinations with up to two par- 



ticles excited. 



we show an example of the Josephson regime, Ui = C/o/2, 
f/oi = (3/4) ?7o, Co = 10^ Ji = 5 Jo, AS = gNUq, with 
iV = 10 atoms and ct = 2 • 10^. In Fig. ^a.) we show 
all the Fock-space amplitudes \ci\^ for all the eigenvec- 
tors, while Fig. [2jb) shows a zoom of these for the lowest 
excited 3A^ + 3{N - 1) -t- 1 states. The first iV 1 are 
eigenstates with no coupling to Fock vectors with atoms 
in the excited level. Then, the following 2N are eigen- 
vectors where the Fock vectors show — 1 atoms in the 
lower level and one atom in the excited level. Finally, the 
following 3(A^ — 1) show two atoms excited and iV — 2 in 
the lower level of eigenenergies (see table |l]). In Fig.[2];c) 
we show the Fock-space amplitudes |cip for the ground 
state, (A^-l-l)th, and (iV-|-2)th excited eigenstates. These 



two figures show the HO behavior discussed in Sec. IV B 
In Fig. [2jd) we represent all the eigenvalues while in the 
inset in this figure we represent a zoom for the lowest 
excited N +1 eigenvectors, where the linear behavior de- 
scribed by Eq. (37) is shown. Every step in the values of 
the eigenvalues in Fig. [2l 



occurs when the correspond- 
ing eigenstate is the superposition of Fock vectors with 
one more atom occupying the second level of energies. 
The increase in the slope in every step is due to the fact 
that Ji> Jq. Indeed, for the last A'' -I- 1 eigenvectors, for 
which all atoms are excited to the higher level, the slope 
is Ji. 

If we reduce Co to one, we cannot assume that the pa- 
rameters correspond clearly to the Fock or to the Joseph- 
son regime, as represented in Fig. [Sja). The Fock-space 
amplitudes \ci\^ of the iV + 1 lowest excited eigenvectors 
show different behavior for the ground state and for the 
{N + l)th excited state, as represented in Fig. [3](b) and 
(c). The ground state is still a HO eigenstate, while the 
(A^-l-l)th excitation is a NOON-like state, as correspond- 
ing to the Fock regime. Fig. [3](d) shows a different be- 
havior than in the previous case. Particularly, as shown 
in the inset in Fig. [sjd), the eigenvalues behave linearly 
for the five lowest excited eigenvectors. Then, the next 
six appear in quasidegenerate pairs with small splittings, 
as corresponds for the Fock regime. Notice that, since 
Ji > Jo, the less excited -I- 1 eigenvectors behave as 
HO eigenstates. Then in this mixed regime, states with 
atoms occupying only the lower level can behave as in 
the Fock regime, while states with atoms occupying the 
excited level behave as in the Josephson one. 

If we decrease C further, to C = 10"^, we enter 



clearly in the Fock regime, as represented in Fig. [4] 
Now, the eigenstates are MS states, as shown in 
Fig. Qb) and (c), while the ground state is sharply 
peaked in |A^/2, iV/2)|0, 0). Also, Fig. gd) shows 
that the MS states appear in quasidegenerate symmet- 
ric/antisymmetric pairs, with small splittings in the 
eigenvalues (see also inset in this figure). Notice that, 
when atoms are excited, the presence of the term 
2C/oi X]j i^t' "tT-jifi-ji' breaks the degeneracy in the Je = 
case between the vectors \Nq — i',i/)\Ni — v\v') and 
|A^o ~ v,v)\v\Ni — v') . For small Jf, this in turn makes 
the MS states with iVi ^ appear also in pairs, and not 
as the superposition of four Fock vectors. These pairs 
can be observed in Fig. |4]^b). 

Finally, the last regime is characterized by the pres- 
ence of shadows of the MS states. Using perturbation 
theory we can show that the MS states Irii^o, "-fl,o)|0, 0)± 
|'^i?,o, '^L,o)|0, 0) couple to MS states showing two ex- 
cited atoms - 2, nfl^,o)|2, 0) ± \nRfi,nLfi - 2)|0, 2) 
and |nL^o,"fl,o - 2)|0,2) ± Iti^^o - 2, nL,o)|2, 0) with co- 
efficients 



^2n^^o(n^^o - 1) 
C/o(6-4ni,o) + 2?7i +2A£; ' 

y^2nfl,o(nfl^o - 1) 
-[/o(6-4nfl,o) + 2C/i + 2A£: ' 

For the most excited MS state, n^^o ~ N, and 

^2N{N-1) 



CL = Uqi 

cr = Uqi 



cn = Uqi 



Uoi6-m) + 2Ui+2AE 



(44) 



(45) 



If this coupling is not negligible, shadows of the MS states 
are coupled, and the upper level cannot be neglected, as 
discussed below. 



C. Crossings of the Eigenvalues 

Let us show, for both regimes, bounds on the one- 
and two-level approximations other than criterion ( 45 1 



With our choice of indexing states, the one-level approx- 
imation corresponds to truncating the size of the Hilbert 
space to -I- 1. With this definition of the one- level ap- 
proximation, the bounds we present below will describe 
the regime in which this truncation is valid. 

Let us consider first the crossings in the Josephson 
regime. The energy levels are coupled by the interac- 
tion energy Uqi. Energy levels only become completely 
decoupled when Uqi = 0. In this regime Uqi is small, and 
then the coupling between levels is weak. Criterion ( 45 ) 



gives a good estimate of this coupling. Let us show that, 
although the energy levels are weakly coupled, eigenvalue 
crossings are induced by the presence of the excited level. 

We consider first Uq = 0, Ui = and Uqi = 0. We find 
that the maximum of the eigenvalues for the first A^ + 1 
eigenstates with no occupation of the excited level coin- 
cides with the minimum of the eigenvalues of the states 



10 



0.05 



0.1 



(a) 



0.15 



0.25 



(a) 



0.2 



0.4 



0.6 



0.8 



250 



200 



150 



100 



liil, 



I III 




N'=10 
N =9 

N'=7 



N'=4 



50 100 

(b) 



200 250 k 

(c) 



i:>i' 1^1 I : I i 

30" 'Tj^^JilT^i 

20 

10 20 30 40 50 k 





— k=0 
--k=10 
— k=11 



20 



30 i 




FIG. 5. (Color online) First crossing in the Josephson regime. 
Panel layout as in Fig. [2] The first crossing in the Josephson 
regime is induced by the reduction of AE or increasing Ji — Jq. 
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FIG. 6. (Color online) First crossing in the Fock regime. 
Panel layout as in Fig. |2] The first crossing in this regime 
is induced by the reduction of AE or the increase in Uq. 



with one atom in the excited level if 

AE 

" Jo i^N - 1) + Ji 



= 1 



(46) 



Equation (46), determines the first eigenvalue crossing in 
this regime. For XJos > 1 no crossing occurs. In Fig. [5] we 



show an example where all the parameters are the same 
as in Fig. [2j that is Co — 10^ , but with a smaller gap be- 
tween levels, since a — 235. As shown in Fig.[5](b) and (c) 
the (7V-|-l)th eigenvector shows occupation of the excited 
level. Notice that for small interactions, the presence of 
the excited level can induce this crossing. Fig.jsjjd) shows 
that the steps in the eigenvalues due to the excited atoms 



11 



observed in Fig. [2|^d) are now absent, since the excited 
atoms have now a much smaller energy. On the other 
hand, the inset in Fig. [2jd) shows that the behavior of 
the eigenvalues is still linear for all the eigenstates which 
are formed by superpositions of Fock vectors with the 
same number of excited atoms. For 



(a) 



XJo 



AE 

Ji — Jn 



= 1 . 



(47) 



the first crossing involving the ground state occurs, i.e., 
the ground state is a state with non-zero occupation of 
the excited level if xjos.gs < 1- This is, indeed, a limiting 
criterion for the validity of the model. 

Analogously, in the Fock regime, the first eigenvalue 
crossings occurs when the condition 



XFock — 



AE 



Uo{Ny2 + N-2)~ 2Uoi{N/2 - 1) 



(48) 



is met. This condition is obtained equating the maxi- 
mum eigenvalue for states with no occupation of the ex- 
cited level, to the minimum eigenvalue for states with one 
particle in the excited level. For XFock > 1 no crossing 
takes place. In Fig. [6] we show an example for which all 
the parameters are the same as in Fig. [4j corresponding 
to the Fock regime, but now cr = 5.1. Here, the crossing 
appeared in the Fock regime, which can be due to the 
reduction of the gap between levels or to the increasing 
of the interactions. Figure [5jd) and its inset show that 
the eigenvectors still appear as quasidegenerate symmet- 



ric/antisymmetric MS states, as discussed in Sec. IV C 

If we increase further the interactions, or reduce the 
energy gap, shadows of the MS states appear, as dis- 



cussed at the end of Sec. |VB| For example, if we reduce 
the AE by taking a = 2.5, the coupling between the MS 
states with no excited atoms and MS states with excited 
atoms is not negligible, as shown in Fig. ^h) and (c). 

Finally, the first crossing involving the ground state 
occurs when the condition 



XFock, 1 



Uo{N 



AE 

-2Uoi{N/2- 



(49) 



is satisfied. For XFock,gs > 1 the ground state shows non- 
zero occupation of the excited level. This criterion should 
be also considered as a limit for the validity of the model. 
As discussed in Sec. |V A[ when this limit is reached one 
must consider more energy levels or an MCTDH p5H5T| 
or other many-body method f47! . 

In Fig. [8) the first iV + l + 2(iV-l) eigenvalues are 
represented when <^ is reduced, for a = 100. The eigen- 
vectors behave as HO-like states for certain region, then 
as MS states, and finally, the most excited eigenstate 
with occupation only of the lower level crosses the least 
excited state with one atom excited when criterion ( 48 ) 
is reached. If ( is increased further, shadows of the MS 
states appear, and eventually, the limit stated by crite- 
rion (49) is reached. 



Therefore, the bounds to the one level approximation 
are given by criteria ( 46 1 and ( 48 1 . We are interested 
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FIG. 7. {ColoT online) Shadows of the MS states. Panellayout 
as in Fig. [2] If Uo is increased, the MS states are coupled 
to Fock vectors with occupation of the excited level. The 
eigenvalues still behave as in the Fock regime. 



in describing cat-like states which are typically excited 
eigenstates. Therefore, a characterization of eigenvalue 
crossings of energies other than the ground state are also 
relevant. These crossings appear when criteria (47) and 
(49 1 are satisfied. 
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FIG. 8. (Color online) Crossings as ( is decreased. Shown are 
the eigenvalues of the A''+l eigenvectors with no occupation of 
the excited level (solid line) and the 2{N— 1) eigenvalues of the 
eigenvectors showing only one excited atom. The ground state 
is non-degenerate, while the rest appear in quasidegenerate 
pairs. The dash-dotted line corresponds to the critical value 



of at which the first crossing occurs, as calculated with ( 48 1 



As (^0 is decreased, the eigenstates change from the Josephson 
regime, where they have are HO states, to the Fock regime, 
where they are MS states. If is decreased further, crossings 
appear, and eventually, shadows of the MS states. 



VI. THREE-DIMENSIONAL DOUBLE WELL 



The three dimensional (3D) double well was exten- 
sively studied in [21]. In this case w ^ uj±. The field 
operators in Hamiltonian (|6| can be expanded in a fixed 
well localized single-particle basis obtained from the delo- 
calized eigenfunctions of the single particle Hamiltonian 

Hsp — 



2M 



V + V^(x). These eigenfunctions are dis- 
tortions of the eigenfunctions of the harmonic oscillator 
potential V^(x) ~ ^w^x^, given by 



(50) 



with n e {0,1,2,...}, i G {n, n - 2, n - 4, . . . , 4iin}, 
and m G {-£, -£+1, £ - I, £}■ Here Rni{r) is the 
radial part of the wavefunction, Yim{0^ </>) are the familiar 
spherical harmonics, and £niin is when n is even and 1 
when n is odd. Then, n is the single-particle energy level, 
£ is the orbital angular momentum in 3D, and m is its 
^-projection. For the two level approximation n = £, 
and therefore the index n is superfluous and hereafter 
suppressed. These functions can be explicitly written in 



spherical coordinates for £ = 0, 1 as 

^oo(r,0,^) = a,-^/\-3/4^-'-V2.L, 

^io(r, 9, ^) = a-„'/'N/27r-3/4(r/aho)e-'-'/2"'" cos(0) , 

^i±i(r, 6, ^) - a-l'\-^'Hrla^o)e-'"'^''-° sin(0)e±''^ . 

The energy of an atom associated with the wavefunction 
V'£„i(x- Xj) is 



En ~ huj{n + 3/2) 



(51) 



In general, the spherical harmonics cannot be accurately 
used as a basis. In particular, they give poor approx- 
imations of the overlap between functions localized in 
different wells. Then, the localized functions at well j 
can be obtained numerically as 



V'joo(x) = (j>{xf(l){y)°tpjo{z) 

V^,io(x) =0(a;)V(y)Vji(^) 
1 

71 

± ^<|>'>{x)<j>\y)^P,o{z)] , 



V'j,i±i(r,0,^)(x) = — [(/.i(a;)0O(y)V',o(^) 



(52) 
(53) 

(54) 
(55) 



where 4>^{x) and (t>^{y) are two lowest excited eigenfunc- 
tions of the harmonic oscillator in the x and y directions, 
and '4>je{z) are the on- well localized eigenfunctions of the 
double well given in (fTol) and (fTTl), respectively. The 



four modes localized in one of the wells are represented 
schematically in Fig. [ijb) . 

Then, we can expand the field operator in ([T]) in terms 
of this eight mode basis as 



*W = ^j^mV'^m(x - Xj) , 



(56) 



where xi = — Xmin and X2 = Xmin are the minima of the 
left and right wells. The operators b^j^^ and bjg,n satisfy 
the usual bosonic annihilation and creation commutation 
relations, 



(57) 



By using this eight-mode expansion of the field opera- 
tor, the 3D Hamiltonian was obtained: 



H — ^ Him + -ffint • (58) 
The first term is a sum over LMGHs for each level, similar 



to the one discussed in Sec. HI [See Eq. ([13|)]: 



Him — ^ njem{njem ^ 1) ^ Jim ^'^jlm^j'lm 

+ Et'Y fljlm ■ (59) 
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In Eq. ( 59 ) , the term Ei J2j njem accounts for the energy 
of the atoms at level £ and z component of angular mo- 



coefRcients in Eq. ( 59 1 are defined as 



mentum m, with hjem 
The term 

Hi, 



^jirrfijim number operator. 



(60) 



accounts for processes among atoms in different levels 
and with different z component of the angular momen- 
tum. 



There are three relevant processes in Eq. (60). The 
first of these is given by 



A°„ter {^01 feoo)'(^io)' + h-C 



(61) 



This process describes the excitation of two atoms from 
the ground state to an orbital with m = 0, or conversely, 
their decay from an excited state with m = to the 
ground state. We name this process zero-vorticity inter- 
level hopping. The second is 



A'nter =E {^01 (^loo) ^ ^11^,1,-1 + h-C 
j 

+ 4:UqI {hjoo {hjii + nji__i))| . 



(62) 



This process is similar to the first one and we call it 
vortex-antivortex interlevel hopping. It permits an atom 
to change its level and also its z component of its angular 
momentum, m. Through this process two atoms in the 
ground state can be excited, one with m = 1 (a vortex) 
and the other with m = — 1 (an anti- vortex) . Conversely, 
one atom with m = +1 and another with m = —1 can 
decay to the ground state. Finally, the third process is 



iJintra =E {^^l feio) '^,11^,1.-1 + h-C 



(63) 



+ 2U°l {hjio (fijii + n^i _i) + 2C/ii (rijiinji _i)| . 

Through the third process, vortex-antivortex intralevel 
hopping, the atoms can only change their angular prop- 
erties, but not their energy level. Then, two excited 
atoms with m = can generate a pair of atoms, one 
with m = 1 and the other with m = — 1 or vice versa. 
These three processes are represented schematically in 
Fig. [TJb). Again, in obtaining this Hamiltonian, the off- 
site interaction coefficients were neglected, as discussed 
in Sec. IIIII for the ID case. 

Notice that the 3D interaction and tunneling coeffi- 
cients that appear in the LMGH part of this Hamiltonian, 
Eq. (59), have to be evaluated independently for atoms 
in the same level and with the same m. The tunneling 



Jlm = - j d^^lPlmi^ - Xmin) 
X V'^m(x -I- Xmin) . 

The interaction coefficients are 



2M 



-f F(x) 



ttI m 

U — 



(fx\^Je„i{:>i-)\^\t/Ji',n'{:>i-W 



(64) 



(65) 



Notice that the coefficients in the coupling part of the 
Hamiltonian, Eqs. (61)-(63), are related to the interac- 



tion coefficients, and not to tunneling coefficients. 

The presence of many new coefficients makes neces- 
sary a wider characterization of regimes, further away 
from the ID Josephson, intermediate, and Fock regimes 
(where the eigenvectors and eigenvalues behave as in 
Figs. [2j[3] and [4]). This was accomplished in [24j, where 
criteria for the crossings and limits of the model were 
also given, and represented numerically for the partic- 
ular case of the Duffing double well potential V{z) = 
Vb(-8zVa^ + 162^/a4 + l). A Josephson and Fock regime 
can be described, for which the eigenvectors can be char- 
acterized as HO-like or MS states, respectively. These 
MS states can show occupation of the excited level, re- 
sulting in an MS of angular degrees of freedom. The 
presence of these orbitals are the fundamental difference 
from the ID case. Moreover, the three coupling pro- 
cesses described above give rise to interesting dynamical 
phenomena, like vortex tunneling and vortex/antivortex 
creation/annihilation along with the tunneling of atoms 
with non-zero m. These new dynamical phenomena are 
the subject of our upcoming work |48) . 



VII. CONCLUSIONS 

Ultracold bosons in double well potentials are a simple 
system to study a great variety of physical phenomena. 
There are two main processes: the atoms can interact in 
pairs on-well with energy Uq or tunnel to the other well 
with energy Jq. When the interactions dominate over the 
tunneling energies in the system, the spectra of the eigen- 
vectors are characterized by the presence of MS states. 
But also in this case it is necessary to consider the pos- 
sibility that the atoms populate an excited level. Then, 
other energies are relevant, including the interaction en- 
ergy of the atoms in the excited level Ui, the tunneling 
energy in the excited level Ji , and the single-particle en- 
ergy gap between levels, AE. 

In this Article, we used a two level approach to de- 
scribe the possible physical scenarios in one- and three- 
dimensional double wells. For the ID case, we clearly 
identified the Josephson regime, for which quantum tun- 
neling was experimentally demonstrated. We character- 
ized the eigenvectors and eigenvalues in the noninter- 
acting regime, and we showed through direct diagonal- 
ization of the two level LMGH that this behavior can 
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be extended to the Josephson regime. Our interest was 
mainly in the occurrence of MS states of atoms locaUzed 
in either one or the other well, which appear for big- 
ger interactions, in the Fock regime. We characterized 
the appearance of crossings in the spectra, which serves 
as a limiting criterion for the validity of the one-level 
or LMGH approximation. We obtained also the limits 
of the two-level approximation. In particular, we found 
when the interactions are sufficiently large that coupling 
to states with non-zero occupation of the excited level 
is relevant. In this case, the interaction coefficient can 
overcome the energy gap between levels, and then the 
states show non-zero occupation of the excited level. We 
described numerically the transition from the Joseph- 
son regime to the Fock regime. We introduced a new 
Josephson-Fock mixed regime, in which the most excited 
eigenvectors with no occupation of the excited level be- 
have as MS states while the less excited ones are HO-like 
states. Moreover, since Ji ^ Jo, when the eigenvectors 
involve only atoms in the excited level, they behave as 
HO-like states. 

MS states are, unfortunately, highly excited states and 
fragile against decoherence processes [371 SH] ■ Thus they 
are difficult to observe experimentally. Quantum super- 
positions of matter waves have been observed for few par- 
ticles, like electrons, but remains a challenging problem 
for larger objects; experiments with Cgo molecules are 
in the lead at present in such efforts (see [5D1 EI] and 
references therein) , but ultracold bosons have the poten- 
tial to go to hundreds or thousands of particles in an 



MS state. Therefore, there are many theoretical propos- 
als for realizing them in a BEG experiment [151 EH 1521 - 
and this remains an appealing research topic with 
deep physical implications. The two-level scenario intro- 
duces the possibility to study other initial states with 
non-zero occupation of the excited level, which could be 
the key to realize MS states experimentally, as discussed 
in [34| . Also, the two-level approach in three dimensions 
permits one to consider angular degrees of freedom in 
the problem [24]. The angular degrees of freedom are a 
very interesting topic in systems of ultracold atoms in 
optical lattices [55M5] , where the excitation of part of 
the population to ap level was demonstrated experimen- 
tally [69] EQj . Here, we summarized the derivation of the 
Hamiltonian and the MS states. Then, the study of the 
dynamics of ultracold atoms in three dimensional double 
wells where orbital degrees of freedom play a relevant role 
is an interesting topic for future research. 
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